subroutine calcwmbej(iijka, iabij, &
          t1, t2, r1, r2, &
          wmbej, icore) 
use mod_iop 
implicit none
real*8, intent(inout) :: t1(*), t2(*), r1(*), r2(*) 
real*8, intent(inout) :: iijka(*), iabij(*)
real*8, intent(inout) :: wmbej(*)
real*8, intent(inout) :: icore(*)
!
!============================
if(ldebug) then 
   call printdnrm2('wmbej_0', wmbej, 4)
endif
call wmbej_t1bn_imnje(wmbej, iijka,  &
                 t1, icore, 1.d0) 
if(ldebug) then 
   call printdnrm2('wmbej_1', wmbej, 4)
endif
call wmbej_iefmn_t1fj_t1bn(wmbej, iabij,&
              t1, icore, -1.d0)
if(ldebug) then 
call printdnrm2('wmbej_2', wmbej, 4)
endif
!============================
call wmbej_iefmn_t2bfjn(wmbej, iabij,&
              t2, icore, 0.5d0)
if(ldebug) then 
call printdnrm2('t2', t2, 3)
call printdnrm2('iabij', iabij, 3)
call printdnrm2('wmbej_3', wmbej, 4)
endif
!============================
call r2abij_t2aeim_wmbej(r2, t2, wmbej, icore, 1.d0) 
!============================
!
return
end 
